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ABSTRACT 

An  investigation  was  made  into  the  rate  of  equilibration 
of  the  stable  isotopes  of  oxygen  in  the  water- sulphate  system 
in  the  ocean.   Two  cores  from  the  Carmel  Submarine  Canyon  off 
the  California  coast  were  analyzed.   The  absolute  ratios  of 
the  oxygen  in  the  interstitial  water  and  dissolved  sulphates 
were  determined  using  a  mass  spectrometer.   The  various  bio- 
logical and  physical  processes  which  can  alter  the  isotopic 
ratios  were  discussed  and  a  simple  mathematical  model  was 
proposed.   The  mathematical  model  retained  terms  which  were 
considered  significant  in  the  anaerobic  zone  of  a  sediment. 
After  simple  manipulations,  a  computational  model  was  derived 
and  a  time  constant  for  the  oxygen  ratio  of  the  water-sul- 
phate system  was  determined  as  -0.0150±  0.0029  (standard 
deviations)  year"1.   From  this,  the  ta  was  determined  to  be 
46.2  years. 
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I.   INTRODUCTION 

All  elements  in  nature  progress  from  states  of  lower 
entropy  to  states  of  higher  entropy.   Along  the  way  to  its 
final  resting  state,  an  element  may  be  acted  upon  by  many 
energy  sources  and  be  cycled  through  a  nearly  closed  loop 
many  times.   At  each  step  of  these  loops,  the  element  is 
acted  upon  while  being  recombined  into  its  next  form  in  the 
loop.   Oxygen  exists  naturally  as  three  stable  isotopes  of 
mass  weight  16,  17  and  18.   Of  these,  017  is  found  only  in 
trace  amounts  so  will  not  be  considered  further.   As  the 
stable  isotopes  of  oxygen  move  from  one  chemical  combination 
to  the  next,  the  ratio  of  one  mass  number  to  the  other  may 
vary  accordingly.   Some  combinations,  notably  C02 ,  rapidly 
equilibrate  with  a  new  environment.   Others,  like  N032", 
P033"  and  S0.+  2",  take  a  relatively  long  time  to  exchange 
their  oxygens  and  exhibit  the  oxygen  isotopic  ratio  of  their 
surroundings.   The  purpose  of  this  paper  is  to  investigate 
the  rate  at  which  sulphates,  derived  from  a  terrestrial 
source,  come  to  equilibrium  with  the  oxygen  in  ocean  waters. 
If  such  a  rate  could  be  determined  for  the  oceans,  it  would 
provide  a  way  to  date  water  masses  without  resorting  to 
radioactive  nuclides.   The  equation  for  this  exchange  is 

H2016  +  S0^6018    t     H2018  +  SO,,16 
The  ratios  of  the  018  to  016  are  defined  as  the  oxygen  ratio 


A.    TECHNIQUE 

Originally  it  was  hoped  that  it  would  be  possible  to 
identify  the  water  mass  emanating  from  the  Carmel  River, 
California,  follow  it  as  it  admixed  with  the  ocean  and  then 
determine  a  rate  of  equilibration  by  knowing  the  age  of  the 
water  mass.   However,  this  method  presented  insurmountable 
difficulties  in  following  a  specific  flood  of  the  Carmel 
River  any  further  than  the  first  headland  so  it  was  decided 
that  the  sediments  of  the  Carmel  Submarine  Canyon,  which  be- 
gins in  Carmel  Bay  (Figure  1) ,  offered  a  handy  source  of  a 
time  series  of  sufficiently  long  duration  so  that  the  time 
rate  of  equilibration  of  the  water- sulphate  system  could  be 
determined  from  the  interstitial  water  and  sulphate  of  a 
standard  core.   While  there  were  many  problems  associated 
with  the  study  of  oxygen  isotope  ratios  in  the  interstitial 
water-sulphate  system  of  a  core,  it  seemed  a  feasible  choice 
Carmel  Canyon  was  chosen  for  study  as  it  is  relatively  un- 
explored compared  to  the  more  familiar  Monterey  Submarine 
Canyon  in  Monterey  Bay,  just  North  of  Carmel  Bay,  and  the 
discharge  of  the  Carmel  River  into  it  is  relatively  straight 
forward.   The  mouth  of  the  Salinas  River  stagnates  most  of 
the  year  and  thus  greatly  complicates  the  history  of  the 
sulphates  it  discharges  into  the  head  of  the  Monterey  Canyon 
Therefore  two  gravity  cores  were  taken  at  the  head  of  the 
Carmel  Canyon  (Figure  1) .   The  interstitial  water  of  these 
cores  was  analyzed  for  the  oxygen  ratios  (the  ratio  of  the 
018  to  the  016)  of  the  water  and  the  dissolved  sulphates. 
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Figure  1.   Locations  of  two  cores  taken  in  the  Carmel  Canyon 
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These  are  reported  herein  as  R„  n   and  R0„   respectively.   In 
the  literature  absolute  ratios  of  this  type  are  not  usually 
reported;  rather,  the  ratio  of  the  sample  is  compared  to  a 
standard  of  known  isotopic  composition.   This  is  done  be- 
cause modern  mass  spectrometers  used  to  determine  isotopic 
abundances  are  of  the  double  collector  type  which  compare 
voltages  from  both  collectors  (unknown  and  standard  samples) 
simultaneously.   The  voltages  are  thus  directly  convertible 
to  relative  ratios.   This  technique  gives  results  accurate 
to  four  and  five  significant  figures.   The  mass  spectrometer 
used  for  the  analyses  in  this  investigation  was  of  the  single 
collector  type  and  therefore  isotopic  abundance  data  is  re- 
ported as  the  absolute  ratios  accurate  to  three  significant 
figures . 

B.    HISTORICAL 

1 .   Laboratory  Studies 

Few  workers  have  addressed  themselves  to  the  problem 
of  the  rate  of  exchange  of  oxygen  isotopes  in  sulphates. 
Teis  (1956),  Hoering  and  Kennedy  (1956)  and  Lloyd  (1967) 
tried  to  estimate  and  explain  the  rate  of  exchange  of  oxygen 
in  sulphates.   Teis  estimated  the  rate  of  exchange  to  be  such 
that  complete  equilibration  occurred  in  about  350  years  (5 
half  lives) .   He  did  not  indicate  the  pH  of  his  solutions 
but  Lloyd  presumed  it  to  be  low  (about  3.2).   Lloyd  and  Teis 
both  let  enriched  sulphate  solutions  stand  at  elevated  tem- 
peratures.  Lloyd  worked  with  higher  pH ' s  and  let  his  samples 
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stand  for  up  to  2  years.   His  data  yielded  an  extrapolated 
time  of  equilibration  of  about  30,000  years  at  4°C  and  a  pH 
of  8.2.   Hoering  and  Kennedy  studied  the  exchange  rate  in 
strong  sulphuric  acid  solutions  (If  to  16f)  where  sulphuric 
acid  is  dissociated  into  the  proton  and  bisulphate  ion.   They 
attempted  to  measure  the  rate  of  exchange  as  a  function  of 
acid  concentration  (activity) .   They  demonstrated  that  the 
exchange  rate  is  a  function  of  pH  and  further  that  the  dis- 
sociation of  H2SO4  into  H20  +  S03    is  the  rate-determining 
step  of  the  equilibrium.   This  is  a  brief  summary  of  the  lab- 
oratory work  which  has  been  done.   All  of  it  has  been  done 
at  temperatures  much  above  the  mean  temperature  of  the  ocean 
(4°C). 

2.   Field  Work 

Some  field  work  by  Rafter  and  Mizutani  (1967)  has 
been  done  to  determine  del  values  for  oxygen  in  the  dissolved 
sulphates  and  the  water  in  Lake  Vanda,  a  highly  stratified 
saline  lake  in  Antarctica.   Their  results  are  reproduced  as 
Figure  2.   Del  values  are  reported  relative  to  SMOW  (Standard 
Mean  Ocean  Water)  (Craig,  1957)  where: 
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and  the  standard  is  sea  water  whose  del  value  is  zero.   If 
some  dates  were  available  for  the  various  layers  of  the  Lake 
Vanda  samples,  a  time  constant  could  be  calculated  directly. 
It  is  suspected  that  the  stratified  portions  of  the  lake 
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Figure  2.   Oxygen  isotope  ratios  of  Lake  Vanda  water  and 
sulphates  versus  depth.   (After  Rafter  and 
Mizutani ,  1967) . 
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have  remained  non-convective  for  up  to  50,000  years  (Clark 
et  al ,  1969)  .   The  samples  Rafter  and  Mizutani  analyzed  had 
been  collected  five  years  earlier  by  Wilson  and  Wellman 
(1962)  and  one  might  suspect  that  in  that  time  span  some 
equilibration  could  have  occurred.   However,  Longinelli  and 
Craig  (1967),  speaking  of  sea  water  samples  from  the  Pacific, 
could  find  "no  tendency  for  samples  stored  up  to  six  years 
to  change  significantly  in  isotopic  composition."   It  is  pre- 
sumed that  the  deeper  layers  represent  progressively  older 
waters.   The  water  then  becomes  progressively  heavier  in  018 
with  depth. 

Teis  (1954)  determined  the  isotope  effect  for  cer- 
tain crystals  and  their  water  of  crystallization.   He  found, 
as  one  might  expect,  that  the  heavier  O18  was  preferentially 
taken  up  by  the  crystal  lattices  of  gypsum,  barium  chloride 
and  sodium  sulphate.   Bi-sodium  sulphate  showed  little  ef- 
fect.  Overall,  one  would  then  expect  highly  stable  struc- 
tures as  the  stratified  layers  of  Lake  Vanda  to  show  a 
tendency  to  become  lighter  with  age. 

The  sulphate  oxygen  in  Lake  Vanda  also  became  pro- 
gressively heavier  with  age  with  a  linear  slope  approximately 
equal  to  0.2  °/00/ft.   Thus  the  water- sulphate  system  is 
being  progressively  depleted  of  O16  by  some  process  other 
than  that  Teis  investigated.   Biological  and  physical  pro- 
cesses as  well  as  time  exchange  have  probably  operated  on 
this  water.   One  Chlorella- like  alga  was  collected  in  these 
layers  (Wilson  and  Wellman,  1962)  and  it  is  expected  that  a 
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commensal  bacterial  population  exists.   Eventually,  if  the 
sulphates  and  the  ambient  water  are  to  equilibrate,  they 
should  reach  the  same  del  value  or  values  whose  ratio  is 
constant.   This  would  give  a  line  of  slope  equal  to  one  on 
a  graph  like  Figure  3.   Such  a  line  has  been  drawn  on  Figure 
3  as  an  extrapolation.   It  can  be  seen  that  the  rate  of 
change  of  the  sulphate  is  greater  than  that  of  the  water. 

van  Donk  and  Mathieu  (1969)  analyzed  water  samples 
from  the  American  Basin  of  the  Arctic  Ocean.   Below  300  m 
the  waters  of  the  basin  are  apparently  non-convective  and 
the  salinity  and  6018  (H20)  are  nearly  constant  with  depth. 
The  water  approaches  the  isotopic  composition  of  SMOW  at 
great  depth.   Unfortunately  no  analyses  were  performed  on 
the  sulphates. 

Rafter  and  Mizutani  (1967)  also  investigated  the 
oxygen  values  of  a  hot  geothermal  spring  at  Wairakei,  New 
Zealand.   This  data,  reproduced  on  Figure  4,  shows  that  the 
water  had  a  constant  del  value  for  all  samples  run  of 
-  5.0  °/oo  (relative  to  SMOW)  and  the  sulphates  had  nearly 
the  del  value  (relative  to  SMOW)  of  ±  2.0  °/0o-   Rafter  and 
Mizutani  assumed  that  this  represented  isotopic  equilibrium 
between  the  sulphates  and  the  water  and  calculated  the  theo- 
retical temperature  to  be  about  140°C.   The  measured  temper- 
ature of  the  bore  was  around  270°C.   They  attributed  the 
discrepancy  to  uncertainties  in  the  theoretical  isotopic 
equilibrium  constants. 
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Figure  3.   Correlation  of  6018  values  in  sulphates  and  the 
associated  water  of  Lake  Vanda,  Antarctica. 
(After  Rafter  and  Mizutani,  1967). 
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Figure  4.   Correlation  of  6018  values  in  sulphates  and  the 
associated  water  of  a  hot  geothermal  spring  in 
Wairakei,  New  Zealand.   (After  Rafter  and 
Mizutani ,  1967)  . 


II.   EXPERIMENTAL  PROCEDURE 

A.  SAMPLING  PROCEDURE 

Two  gravity  cores,  2.5  inch  diameter  and  about  50  inches 
long,  were  taken  from  a  fairly  level  area  of  the  upper  head 
of  the  Carmel  Canyon  (Figure  1) .   The  cores  were  not  extruded 
but  were  kept  in  their  plastic  liners  until  analyzed.   The 
cores  were  frozen  to  prevent  continuing  bacterial  ac- 
tivity and  interstitial  water  migration.   Prior  to  beirtg 
frozen,  the  residual  water  atop  the  core  in  the  barrel  was 
sampled  and  was  also  frozen.   The  two  cores  were  labeled  A- 2 
and  B-l. 

B.  EXPERIMENTAL  METHOD 

One  core  (A- 2)  was  sawn  lengthwise  and  each  half  was 
analyzed  independently  as  cores  A- 2  and  A- 2b.   Both  cores 
(A- 2  and  B-l)  had  approximately  2  inch  lengths  cut  out  of 
them  every  10  inches  giving  six  samples  per  core  plus  the 
water  taken  from  the  top  of  each  which  was  labeled  top-of- 
core.   The  frozen  sediment  was  thawed  and  diluted  with  ap- 
proximately 250  cc  of  distilled  water  and  then  filtered  and 
centrifuged  to  separate  out  the  solids  from  the  effluent. 
Approximately  10  ml  of  0 . IN  NaOH  and  10  ml  of  0 . 1M  A1C13  was 
added  to  the  effluent  to  bring  down  the  suspended  clays  in 
the  water  and  the  flocculant  and  adsorbed  clays  were  centri- 
fuged off.   Approximately  20  ml  of  0 . IN  HC1  and  10  ml  of 
0 . 1M  BaCl2  were  then  added  to  the  effluent  and  the  precipi- 
tate (BaSO^)  was  digested  for  a  half  hour.   The  use  of  a 
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centrifuge  to  wash  the  precipitate  made  digestion  less  neces- 
sary than  if  filtration  had  been  used  alone.  The  precipitate 
was  then  washed  three  times  with  distilled  water  and  dried  at 
90°C  overnight. 

The  BaSC\  yield  from  this  procedure  varied  from  apparent 
volumes  of  5  ml  in  the  surface  layers  of  the  sediments  to 
about  0.2  ml  at  the  lowest  depths  (50  inches)  in  the  sedi- 
ment.  After  drying,  the  precipitate  was  ready  to  be  reduced. 
The  reduction  was  carried  out  as  described  by  Clayton  and 
Epstein  (1967) .   A  Lepel  RF  furnace  was  used  to  drive  the 
reaction.   Figures  5a  and  5b  lay  out  the  experimental  ap- 
paratus.  The  precipitate  was  mixed  with  10  to  20  times  the 
stoichiometric  amount  of  spectrographic-quality  carbon.  Lamp- 
black was  tried  but  residual  organic  compounds  masked  the 
desired  results.   The  crucible  and  pedestals  were  turned 
from  a  rod  of  spectrographic-quality  graphite.   The  crucible 
was  fitted  with  a  tightly  fitting  cover  of  graphite  to  en- 
sure an  even  internal  heat  flux.   A  Leeds  and  Northrup  hot 
wire  optical  pyrometer  was  used  to  calibrate  the  settings  on 
the  RF  furnace  necessary  to  give  temperatures  of  700°C  and 
1150°C. 

The  mixed  precipitate  and  carbon  was  put  in  the  crucible 
and  the  whole  system  was  outgassed  to  about  4x10" 4  mmHg  by  a 
mechanical  pump  backed  by  an  oil  diffusion  pump.   The  RF  fur- 
nace was  turned  on  to  700  °C  to  outgas  all  volatile  products 
and  gasses  which  the  chamber  and  sample  may  have  picked  up 
in  handling.   After  the  pressure  dropped  below  4x10" k   mmHg 
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Cover 

Plalinum    Heat    Shield 


Reaction     Vessel 
Carbon     Support 

Pedestal 

Glass    Support 


Figure  5a. 


Diagrammatic  sketch  of  barium  sulphate  reduc- 
tion chamber.   The  platinum  shield  is  necessary 
to  prevent  the  glass  reaction  vessel  from  melt- 
ing.  The  platinum  shield  must  form  a  complete 
loop  but  must  not  be  closed  so  it  does  not  it- 
self become  heated  by  the  RF  furnace. 
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(about  45  min.) ,  the  pumps  were  separated  from  the  system 
and  the  heat  increased  to  1150°C.   The  reaction  was  allowed 
to  proceed  for  45  min.   The  carbon  dioxide  and  carbon  mon- 
oxide given  off  was  collected  in  an  ampule  which  was  sealed 
with  a  torch  at  the  end  of  the  reaction  period.   The  final 
pressures  recorded  were  about  5  mmHg.   The  top-of-core  water 
was  handled  in  the  same  fashion  since  it  generally  had  a 
large  amount  of  suspended  sediment. 

The  gas  collected  was  analyzed  on  the  Naval  Postgraduate 
School  mass  spectrometer  CEC  type  21-103.   This  spectrometer 
is  a  single  collector  type  and  gives  the  absolute  mass  abun- 
dance ratios  accurate  to  three  significant  figures.   Since 
it  was  not  possible  to  run  samples  against  a  known  standard, 
all  data  is  reported  as  the  ratio  of  the  mass  peaks  at  44 
plus  45  versus  46.   The  mass  ratios  so  derived  are  not  the 
absolute  ratios  of  oxygen  since  the  carbon  atoms  were  in- 
cluded in  the  ratios  calculated.   However  the  oxygen  in  each 
ratio  can  be  easily  determined. 

RSO,  (true)  =  Wm    '    RS0,  (measure-d) 

Each  ampule  (depth)  was  used  for  three  determinations  of  the 
ratio.   These  are  listed  as,  for  instance,  0-2  inches  a, 
b  and  c  in  Table  3.   Three  determinations  were  needed  to 
average  out  transient  noise  spokes  in  the  zero  line  of  the 
mass  spectrogram. 
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Thus, 

Rc„   (measured)  =  ^ - 


it  it   T  r   it  5 


where  P   is  the  value  of  the  mass  peak  at  mass  number  a. 
a  r 

The  mass  peaks  at  44  and  45  constitute  the  sum  of  the  016, 
017,  C13  and  C12  components  and  the  46  peak  represents 
ci20i60i8i   Craig  (1957)  indicated  that,  in  del  values,  the 
correction  for  the  mass  ratio  (C13/C12)  of  the  carbon  used 
is  around  one  percent.   Since  this  is  on  the  order  of  ac- 
curacy of  the  mass  spectrometric  data  used,  it  was  therefore 
not  included  herein. 

In  order  to  determine  the  mass  ratios  of  oxygen  in  the 
interstitial  waters,  the  frozen  sediment  was  first  carefully 
scraped  to  remove  the  exterior  layers.   Then  the  frozen  sedi- 
ment was  put  in  a  vacuum  vessel  and  kept  frozen  by  liquid 
nitrogen  while  the  vessel  was  outgassed.   When  the  pressure 
was  less  then  3x10" k   mmHg ,  the  system  was  closed  off  and  al- 
lowed to  come  to  room  temperature  overnight.   The  water  vapor 
and  other  gasses  given  off  were  collected  and  analyzed  on 
the  mass  spectrometer  like  the  sulphate  samples.   The  water 
ratios  are  reported  as  a  fraction  of  the  ratio  between  the 
mass  18  peak  plus  the  mass  19  peak  versus  the  mass  20  peak. 
Analogous  to  the  sulphate  ratios,  a  correction  must  be  ac- 
counted for  because  of  the  hydrogen  in  the  samples.   This 
can  be  shown  to  be: 

RH20  (true)  =  Wlf  '  RH20  (measured) 
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where  again, 

P2o 

Ru  n  (measured)  =  ~ — 
n  2  U  ru 

The  mass  peak  at  19  was  more  then  100  times  smaller  than  the 
mass  peak  at  18  and,  therefore,  not  included. 

C.  WATER  CORRECTION 

For  a  single-stage  distillation  process  at  15  to  20°C 
the  water  vapor  has  a  6018  of  -8%0and  6D2  of  -\%0.    This  as- 
sumes that  the  evaporation  and  condensation  of  the  water 
takes  place  under  conditions  of  thermodynamic  equilibrium, 
and  the  D2  and  018  content  in  the  liquid  and  gaseous   phases 
is   controlled  by  the  vapor  pressures  as  in  Rayleigh  distil- 
lation (Rankama,  1963).   Therefore  the  absolute  ratio  of  water 
vapor  should  shift  -1.2xl0~5  per  °/00.   This  is  the  correc- 
tion factor  for  water  at  20°C;  however,  since  the  data  in 
this  investigation  was  accurate  only  to  ±  10x10" 5,  it  was 
ignored. 

D.  CHEMISTRY  OF  SULPHATE  ANALYSIS 

As  the  clays  and  silts  are  laid  down  on  the  Carmel  Can- 
yon floor,  water  adsorbed  on  such  particles  and  carried  down 
with  them  is  collected  in  the  sediments  as  interstitial  water. 
Some  of  this  water  is  carried  down  the  canyon  by  the  predom- 
inately down-canyon  transport  from  the  main  fresh  water 
source  in  the  area  -  the  Carmel  River.   Thus  the  sulphates 
in  the  interstitial  waters  of  the  canyon's  sediments  are 
derived  from  the  surface  waters  of  the  California  coast  and 
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the  Carmel  River.   The  derived  water  samples  were  slightly 
acidified  to  insure  that  the  one  barium  salt  to  come  down 
was  BaSOit  according  to  the  equation: 

Ba2+  +  S042"  -*  BaSCn 

Other  salts,  including  borates,  are  too  soluble  to  be  brought 
out  of  solution.   Thus,  oxygen  in  the  sulphates  was  fixed 
as  the  precipitated  barium  sulphate. 

The  v/ashed  precipitates  were  mixed  with  at  least  ten 
times  the  amount  of  carbon  needed  for  complete  reaction.   The 
reduction  was  carried  out  in  a  carbon  crucible.   The  follow- 
ing equations  describe  the  products  of  this  reduction: 

BaSCK  +  C  ->  BaS  +  CO  +  C02  -  HEAT 
C  +  C02  ->  CO  -  HEAT 

Longinelli  and  Craig  (1967)  and  Rafter  (1957)  have  discussed 
the  temperature  dependence  of  these  reactions.   In  general 
lower  temperatures  (1000°C)  of  reaction  produce  more  carbon 
monoxide  whereas  at  higher  temperatures  (1150°C)  less  carbon 
monoxide  is  produced. 

A  series  of  analyses  were  run  on  pure  carbon  dioxide 
and  on  air  to  determine  the  breakdown  of  the  main  component 
and  possible  contaminants  to  the  analysis.   Each  sample  was 
tested  through  the  range  of  operating  pressures  of  the  ion 
chamber  in  the  mass  spectrometer.   The  results  of  these  tests 
are  depicted  in  Table  1.   The  National  Bureau  of  Standards 
published  spectral  data  on  Argon,  C02 ,  CO  and  02  (American 
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Petroleum  Institute  Research  Project  44,  1960)  are  presented 
in  Table  2. 

Table  1 
Mass  Spectral  Data  for  Air  (at  different  pressures) 


Mass  Number 

14 

16    18 

28 

29 

32 

40 

44 

Intensity 

106.5 

31.4 

1248 

9.0 

259.0 

28.3 

- 

Pressure 

19.89 

microns 

Intensity 

195.3 

57.6   2.9 

2394 

17.9 

490 

56.5 

3.5 

Pressure 

17.95 

microns 

Intensity 

309 

90.6   5.6 

3980 

29.6 

825 

98.4 

- 

Pressure 

17.20 

microns 

Intensity 

361 

35.3   5.8 

4790 

35.6 

987 

112.2 

7.9 

Pressure 

20.83 

microns 

Intensity 

394 

116.7   6.3 

5160 

38.3 

1065 

120.3 

8.4 

Pressure 

19.89 

microns 

Intensity 

387 

115.8   6.8 

5000 

37.6 

1041 

117.3 

8.5 

Pressure 

18.97 

microns 
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Table  2 

Published  Mass  Spectral  Data 

ARGON 

Mass  Number     20     36     40 

Relative 

Intensity    13.0   0.32    100 

CARBON  DIOXIDE  #157 

Mass  Number     12     16      28     29     44 

Relative 

Intensity    6.67   9.40    8.17   0.10   100 

#  1582 

Mass  Number    16     22      28     29     44    45 

Relative 

Intensity    5.57   1.36    6.70   0.14   100   1.12 
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0.39 


#  1596 

Mass  Number     12     16      22     28     44    45     46 

Relative 

Intensity    2.46   6.24    1.78   6.55   100   1.16    0.41 


Mass  Number 

Relative 

Intensity 


CARBON  MONOXIDE 
12     14      16     28     29    30 

4.71   0.75    1.67   100    1.16  0.22 


Mass  Number     14     28 


29 


Relative 

Intensity 


5.18   100     0.73 


Mass  Number 

Relative 

Intensity 


OXYGEN 
16     32     33     34 

5.13   100     0.08   0.41 
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III.   RESULTS 

A.    CORE  B-l 

The  sulphate  ratios  with  depth  and  the  overlying  water 
had  a  range  of  RSQ   of  43.1xl0~'t  ±  0.9  (Table  3).   Plotted 
against  depth  (Figure  6)  they  fell  on  a  complex  line.   There 
was  fractionation  apparent  between  the  overlying  waters  and 
the  sediment  water  interface.   The  sulphates  in  the  top  lay- 
ers of  the  sediment  columns  were  rejecting  the  018  and  re- 
taining 016.   There  was  an  inclusion  of  the  heavier  O18  or 
an  exclusion  of  the  lighter  O16  between  the  surface  layers 
and  the  10  inch  depth.   Then  from  10  inches  to  30  inches  in- 
clusive, the  Rcr.   was  constant  at  43.5  ±  0.2xl0_l*.   At  40 
and  50  inches  there  was  a  slight  lightening  of  the  sulphates 
as  the  Rc0   dropped  to  42.7x10  k.      Therefore,  in  this  region 
some  mechanism(s)   was  (were)  operating  to  include  more  of 
the  light  O16  or  to  reject  more  of  the  heavy  O18  from  the 
interstitial  sulphate  ions. 

The  results  of  the  water  analysis  showed  interesting 
variations  with  depth  which  tended  to  correspond  to  the 
sulphate  analysis.   The  top-of-core  water  (Figure  7  and 
Table  4)  representing  overlying  water,  was  lighter  than  that 
from  the  top  sediment  layers,  thus  water  which  came  down  as 
interstitial  water  in  the  sediments  was  preferentially 
heavier  or  was  acted  upon  immediately  to  fractionate  out  the 
O16.   In  the  first  10  inches  of  the  sediment-,  the  water 
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Table 

3 

] 

^or.      Core 

oU  4 

B-l 

DEPTH  IN  CORE 
(Inches) 

(measured) 

Average 

(true) 

Top 

of 

Core 

a. 
b. 
c. 

41.0 
40.6 
41.9 

41.2 

43.8 

0  -  2 

a. 
b. 
c. 

41.6 
40.1 
37.5 

39.7 

42.2 

10 

a. 
b. 
c. 

40.8 
40.8 
40.7 

40.8 

43.4 

20 

a. 
b. 
c. 

41.2 
40.4 
41.0 

40.9 

43.5 

30 

a. 
b. 
c. 

40.9 
40.7 
40.9 

40.9 

43.5 

40 

a. 
b. 
c. 

40.6 
40.5 
38.8 

40.5 

43.0 

-Disregarded 
due  to  noise 
level 

50  a.   41.6  --Disregarded 

40.2       42.7  due  to  noise 

level 


(All  values  of  R   are  presented  xlO1*) 

A. 


a. 

41.6 

b. 

40.5 

c . 

39.9 

30 


Top  of  Core 
00 


Rso4(*1°4) 


Figure  6.   Absolute  oxygen  ratios  in  sulphates  of  core  B-l 
versus  depth. 
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Figure  7.   Absolute  oxygen  ratios  in  water  of  core  B-l 
versus  depth. 
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Table  4 


R 


H20  Core  B-l 


Depth 

in 

Core 

Rr 
(meas 

(Inc 

;hes 

0 

ured) 

Average 

Top 

a . 

22.4 

of 

b. 

21.6 

22.3 

Core 

c . 

22.8 

0  -  2 

a. 

23.5 

b. 

20.9 

23.8 

c . 

24.1 

10 

a. 

22.8 

b. 

23.1 

23.3 

c. 

24.6 

20 

a. 

24.8 

b. 

25.5 

25.1 

c. 



30 

a. 

26.2 

b. 

26.1 

26.2 

c. 



40 

a. 

22.6 

b. 

22.3 

22.7 

c. 

23.6 

50 

a. 

20.9 

b. 

21.0 

21.4 

c. 

22.4 

R 


(true) 
22.5 


24.1- -Disregarded 
due  to  noise 
level 

23.6 


25.4 


26.5 


23.0 


21.7 


(All  values  of  R   are  presented  xlO1*) 
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appeared  to  remain  essentially  constant  or  slightly  lighter 
at  R,,  r,   =    24. 1x10" h.      Then  in  the  next  20  inches  of  core 

112  U 

from  10  to  30  inches,  the  Rir  n  ratio  increased  almost  lin- 

11  2  U 

early  to  26. 5x10"  "* .   In  this  region  some  combination  of 
factors  was  operating  to  remove  the  lighter  O16  from  the 
interstitial  water  or  add  the  heavier  O18.   Below  30  inches 
the  water  became  progressively  lighter  with  depth  falling 
from  a  ratio  of  26.5xl0-lt  at  30  inches  to  21.7xl0_lt  at  50 
inches.   Therefore,  in  this  region  the  heavier  O18  was  be- 
ing taken  from  the  water  or  the  lighter  O16  added  to  it. 
Figure  8  summarizes  the  relative  change  in  O16  with  depth 
for  the  water- sulphate  system  of  core  B-l. 

B.    CORE  A- 2 

Core  A-2  (Figures  9  and  10,  and  Tables  5  and  6)  was 
the  first  to  be  analyzed  and,  due  to  inexperience  in  handl- 
ing and  the  fact  that  the  mass  spectrometer  had  a  higher 
noise  level  than  when  samples  from  B-l  and  A- 2b  were  run, 
had  a  wider  scatter  in  data  points  then  core  B-l.   This 
greatly  complicated  a  discussion  of  results  from  this  core. 
General  outlines  relative  to  core  B-l  could  be  drawn,  how- 
ever . 

The  sulphate  ratio  determined  from  the  overlying  waters 
for  the  two  cores  was  about  the  same  and  a  value  of  43.5xl0_l* 
was  adopted  as  representative  of  the  sulphates  of  the  over- 
lying waters.   The  surface  layers  of  the  sediment  had  about 
this  ratio  within  the  spread  of  data  values  obtained.   Be- 
tween the  surface  layers  and  10  inches  the  sulphate  became 
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AR  (*104) 
0.0     +2.0 


50-J 


Figure  8.   The  loss  or  gain  of  oxygen  versus  depth  for  the 
water- sulphate  system  in  core  B-l.   AR  = 
[(RSQ   (surface)  -  Rg0   (depth)) +  (&H  0  (surface) 
-  RH  Q  (depth)) ]xlO 


-  k 
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Top  of  Core 
00 
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Figure  9.   Absolute  oxygen  ratios  in  sulphate  of  core  A-2 
versus  depth. 
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R 


Depth  in 
(Inch< 

Core 

2S) 

Top 
of 

Core 

00 


10 


20 


30 


50 


52 


SO 


Table  5 
Core  A-2 


R 


(measured) 


a , 
b, 
c 

a 
b 
c 

a 
b 
c 

a 
b 
c 

a 
b 
c 

a 
b 
c 

a 
b 
c 


40.9 
40.6 
40.6 

40.6 
41.8 
41.2 

41.5 
39.8 
39.0 

42.2 
41.0 
41.1 

41.7 
40.8 
41.2 

44.2 
42.5 
42.3 

40.7 
37.8 
38.8 


Average 
40.7 

41.2 

40.1 

41.4 


41.2 


43.0 


38.8 


R 


(true) 
43.2 

43.7 

42.6 

44.0 

43.8 

46.0 

41.5 


(All  values  of  R   are  presented  xlO1*) 

A, 
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Figure  10.   Absolute  oxygen  ratios  in  water  of  core  A-2 
versus  depth. 


Table  6 


R„  _  Core  A- 2 
H2U 


Depth  in  Core     R 


R 


(Inches)        (measured)    Average      (true) 

30  a.  23.8 

b.  24.8      24.9       25.3 

c.  26.3 

52  a.  24.1 

b.  26.0 

c.  25.5       25.5        25.9 

d.  26.6 


(All  values  of  R   are  presented  xlO1*) 

.A. 
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lighter.   Then  between  10  and  20  inches  the  ratio  increased 
to  about  43. 7x10" k    and  was  constant  to  30  inches.   The  bottom 
sampling  depths  of  the  core  were  not  taken  as  representative 
of  the  core  since  it  appeared  that  the  bottom  of  the  core  had 
some  extra  water  frozen  in  with  it. 

The  interstititial  water  analyses  from  A- 2  were  incon- 
clusive because  of  the  lack  of  data  points. 

C.    CORE  A- 2b 

Core  A- 2b  was  the  last  series  of  analyses  run  and  the 
increased  level  of  experience  in  sample  handling  was  reflect- 
ed in  a  tighter  grouping  of  data  points  (Figures  11  through 
13  and  Tables  7  and  8).   Core  A-2  was  taken  about  a  mile  and 
a  half  up  the  canyon  from  core  B-l.   This  is  reflected  in 
the  analysis  of  the  core  which  showed  a  more  complicated 
picture  than  did  B-l.   The  relatively  constant  integrated 
value  of  total  016  was  interpreted  to  mean  that  the  pro- 
cesses which  act  to  remove  or  add  01G  were  acting  approxi- 
mately in  equilibrium  and  therefore  authigenic  mineral 
formation  was  probably  not  occurring.   Over  the  whole  core, 
the  water-sulphate  system  appeared  to  have  lost  a  little 
O16  relative  to  the  surface  (Figure  13).   Coming  from  a  lo- 
cation nearer  the  coast,  a  deeper  aerobic  layer  was  expected. 
Below  40  inches,  the  system  was  probably  beginning  to  be- 
have more  like  B-l,  with  the  water  ratios  moving  to  the  left 
and  the  sulphate  ratios  beginning  to  respond  to  reach  this 
new  equilibrium. 
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Rso  (-104) 
44     4.6 


50-J 


Figure  11.   Absolute  oxygen  ratios  in  sulphate  of  core  A-2b 
versus  depth. 
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Table  7 


Rcn   Core  A- 2b 

oU  it 


Core  in  Depth    R 


R, 


(Inches)         (measured)    Average      (true) 

0-3  a.  41.1 

b.  40.6       40.6       43.1 

c.  40.3 

9-12  a.  42.0 

b.  41.3        41.5        44.1 

c.  41.1 

19  -  21  a.  41.0 

b.  40.1       40.0       42.5 

c.  39.8 

29-31  a. 42.1  This  ampule 

b.  40.9       41.4       44.0   may  have  had 

c.  41.3  air  contam- 

ination. 

39  -  41  a.  40.6 

b.  40.7       40.8       43.4 

c.  41.2 

49  -  51  a.  40.9 

b.  40.2       40.8       43.4 

c.  41.3 


(All  values  of  R   are  presented  xlO1*) 
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Figure  12.   Absolute  oxygen  ratios  in  water  of  core  A-2b 
versus  depth. 
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Table  8 


Ru  n  Core  A- 2b 
H2U 


Depth  in  Core     R  R 

(Inches)         (measured)    Average      (true) 


2  -  3 


10  -  11 


19  -  20 


31  -  32 


41  -  42 


48  -  49 


a . 

b. 

21. 

,6 

c. 

21. 

4 

a. 

b. 

21, 

,5 

c . 

21, 

,4 

a. 

20, 

,7 

b. 

21, 

,1 

c . 

20, 

.8 

a . 

20, 

,7 

b. 

21, 

,  3 

c. 

20, 

,2 

a . 

21, 

,6 

b. 

21, 

,5 

c . 

20, 
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Figure  13.   The  loss  or  gain  of  oxygen  versus  depth  for  the 

water- sulphate  system  in  core  A-2b.   AR  = 

[(RSQ   (surface)  -  RSQ   (depth))  +  (R^q  (surface) 

-  RH  Q  (depth))  JxlO1* 
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IV.   DISCUSSION  AND  THEORY 

A.    GENERAL 

The  underlying  principle  that  one  uses  in  following  a 
suite  of  isotopes  through  a  chemical  system  is  that  the 
lighter  isotopes  react  more  rapidly  than  the  heavier  when 
bond-breaking  predominates.   Thus,  a  process  which  acts  to 
place  the  oxygen  (in  this  investigation)  in  a  different 
energy  state  will  tend  to  fractionate  the  parent  population 
leaving  018  behind  and  removing  0ie,  or  vice  versa  if  bonds 
are  being  formed. 

1 .   Biological  Isotope  Effect 

Microorganisms  (bacteria)  exert  a  strong  influence 
on  the  geochemical  sulphur  cycle  and  it  is  expected  that  they 
will  likewise  influence  the  geochemical  oxygen  cycle  as  it 
pertains  to  sulphates.   There  are  bacteria  which  oxidize 
sulphur  and  its  inorganic  compounds  both  photosynthetically 
(Thiorhodaceae ,  Chlorobacteriaceae  and  Athiorhodaceae)  and 
non-photosynthetically  (Beggiatoaceae ,  Achromatiaceae  and 
Thiobacillus) .   There  are  those  bacteria  which  reduce  sul- 
phates to  hydrogen  sulphide  (Desulf ovibrio  and  Sporovibrio) 
and  there  are  those  which  reduce  organic  sulphur  compounds 
to  H2S  (various  pathogenic  and  saprophytic  types  like 
Proteus ,  Serratia  and  Clostridium) (Frobisher,  1957). 

Of  the  bacteria  which  form  sulphate  from  reduced 
sulphur,  most  require  free  oxygen  to  live.   An  example  of 
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an  aerobic  Thiobacillus  reaction  would  be 

SNa2S203  +  H20  +  402  ->  5Na2S0^  +  tizSOh    +  4S 

which  is  the  overall  reaction  for  Thiobacillus  thioparus . 
There  is  also  an  important  group  of  sulphate  producing  bac- 
teria which  do  not  need  oxygen  but  can  utilize  nitrates  as 
an  energy  source.   A  common  one  found  in  soil  is  Thiobacillus 
denitr if leans ;  its  overall  reaction  is 

5S  +  6HNO3  +  2H20  ->  5H2S04  +  3N2 

Frobisher  points  out  that  the  thionate  bacteria  can  occur 
where  sulphur  and  its  compounds  occur  as  a  result  of  chem- 
ical action  or  due  to  protein  decomposition  of  the  sulphate 
reducing  bacteria  Desulf ovibrio .   It  is  of  some  importance 
here  to  note  that  the thionate  and  nitrophilic   bacteria  use 
water  in  their  reactions  as  well  as  thionate,  oxygen,  sul- 
phur and  nitrates,  thereby  establishing  oxygen  exchange  with 
the  ubiquitous  solvent. 

The  sulphate  reducing  bacteria  are  divided  into  two 
genera.   Sporovibrio ,  however,  is  thought  to  be  only  a  spore- 
forming  variant  of  Desulfovibrio .   The  fresh  water  specie  D. 
desulfuricans  is  a  strict  anaerobe  and  heterotroph.   Its 
overall  reaction  can  be  expressed  as 

2CH3'CHOH-COONa  +  HsSCU  +  2CH3-COONa  +  H2S  +  H20  +  2C02 

Notice  that  two  end  products  are  carbon  dioxide  and  water. 
Some  sulphate  reducers  are  autotrophic  and  can  directly  re- 
duce inorganic  sulphates  by  means  of  the  enzyme  hydrogenase 
4H2  +  H2S0^enzymel  H2S  +  H20 

4  7. 


Again,  water  is  an  end  product.   The  photosynthetic  reactions 
will  not  be  considered  due  to  the  location  of  the  samples. 

There  have  been  a  considerable  number  of  laboratory- 
studies  done  on  the  fractionation  of  sulphur  in  the  biologi- 
cal reactions  discussed  above.   Nakai  and  Jensen  (1964)  in- 
vestigated the  biological  isotope  fractionation  on  sulphur 
in  sulphates  and  sulphides.   They  showed  that  the  isotope 
fractionation  of  sulphur  in  sulphates  was  constant  at  1.020 
when  the  supply  of  sulphate  was  limited  in  a  closed  system. 
Thus,  regardless  of  the  instantaneous  concentration  of  sul- 
phate, S32  will  react  2%    faster  than  will  S34.   They  also 
showed  that  the  reaction  is  unidirectional  and  comparable  to 
that  expected  from  an  inorganic  reduction  of  sulphates. 
Harrison  and  Thode  (1958)  postulated  that  the  bacterial  re- 
duction of  sulphates  caused  fractionation  of  the  sulphur  in 
two  possible  ways,  namely  in  the  reduction  of  SCU2~  to  S032" 
there  is  first  a  combining  step  of  the  sulphate  to  an  enzyme 
and  then  secondly  the  enzyme- sulphate  bond  is  broken  with 
the  enzyme  cleaving  off  an  oxygen  from  the  sulphate;  i.e., 

STEP  I 

..  Cl 

SOt,    +  enzyme   ■•>   SOi^enzyme 

c2 

STEP  II 

c3     , 
enzyme*  SO4   -»■  S03    +  Oenzyme 

They  said  that  step  I  left  behind  S32  in  the  free  sulphate 
ions  and  preferentially  took  up  S31*.  Then  in  step  II,  the 
oxygen  cleavage  step,  the  sulphate  complex  with  S32  was 
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preferentially  broken  to  S032~.   This  postulate  is  in  accord 
with  the  general  principle  stated  originally.   Harrison  and 
Thode  showed  at  very  low  sulphate  concentrations  (1x10  5moles/ 
1)  that  step  I  became  unidirectional  in  the  Cj  direction  and 
dominant,  while  at  higher  concentrations  step  II  was  dominant. 

Biological  reduction  of  sulphates  have,  in  general, 
been  shown  to  fractionate  sulphur  (when  there  is  a  fractiona- 
tion) such  that  the  residual  sulphate  is  enriched  in  S 3  ** 
(Yeremenko  and  Mekhtiyeva,  1961;  Kaplan  and  Rittenberg,  1964; 
Jones  and  Starkey,  1957;  Kaplan,  Emery  and  Rittenberg,  1963; 
Nakai  and  Jensen,  1964;  Mekhtieva  and  Kondrat'eva,  1966; 
Thode,  Kleerekoper  and  McElcheran,  1951)  except  at  very  low 
sulphate  concentrations.   When,  conversely,  sulphate  is  being 
biologically  synthesized  it  will  be  lighter  in  S3lt  than  the 
source  of  sulphur  (Nakai  and  Jensen,  1964). 
2 .   Sulphate  Distribution  in  Sediments 

A  simplified  model  of  the  distribution  of  sulphates 
in  recent  sediments  was  proposed  by  Berner  (1964).   Berner's 
model  assumed  that  the  local  time  rate  of  change  of  concen- 
tration of  sulphate  in  the  interstitial  water  of  sediments 
was  the  sum  of  three  terms: 

a.  Ionic  diffusion 

b.  Deposition  plus  compaction 

c.  Bacterial  (anaerobic)  sulphate  reduction. 

His  solution  was 

-(kx) 
C  =  (Cq  -  Coo)  e  w    +  Coo 
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where , 

c   =  concentration  of  sulphates, 

w   =  rate  of  deposition  less  rate  of  compaction  and 

is  equal  to  dx/dt, 
k   =  carbon  oxidation  rate  constant  for  sulphate 

reducing  bacteria, 

x   =  depth  (in  the  core  in  cms) , 

C   =  initial  concentration  (at  x  =  0) , 
o 

Coo  =  asymptotic  concentration  (at  x  =  °°)  . 

Berner  fitted  his  solution  to  sulphate  concentrations  in  one 
of  the  cores  of  Kaplan  e_t  a_l  (1963)  .   This  is  reproduced  as 
Figure  14.   The  equation  for  the  fitted  curve  is 

C  =  19  e("-015x)  +  8 

Berner  fudged  a  bit  here  since  values  of  sulphate  concentra- 
tion increase  after  depths  greater  than  320  cm.   Kaplan  ejt 
al  (1963)  suggested  that  the  asymptotic  value  for  sulphate 
concentration  at  depth  could  be  an  artifact  of  analytical 
technique  or  could  be  a  measure  of  the  exhaustion  of  utiliz- 
able  organic  carbon.   Nakai  and  Jensen  (1964)  showed  that 
biological  oxidation  of  sulphur  and  pyrite  took  at  least  10 
days.   Thus  Kaplan  et  a_l  were  probably  wrong  to  suspect  mani- 
pulation as  a  source  of  error.   Two  of  the  three  cores  ana- 
lyzed by  Kaplan  e_t  a_l  showed  a  virtual  monotonic  decrease  in 
sulphate  concentration  with  depth.   One  (Santa  Catalina  - 
6248)  showed  a  zone  of  increasing  sulphate  concentration  with 
depth  in  the  first  meter  of  the  core,  then  a  steady  decrease 
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Figure  14.   Sulphate  concentration  versus  depth  data  for 

core  7131  from  Santa  Barbara  Basin  (after  Kaplan 
et  al,  1963) . 
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in  concentration  until  values  were  comparable  to  the  other 
cores  (Figure  15) .   The  Kaplan  group  felt  that  the  surface 
of  the  sediments  is  an  area  where  aerobic  bacteria  are  ac- 
tively fractionating  sulphur  compounds,  giving  the  sediment 
column  a  population  of  combined  sulphur  whose  isotopic  com- 
position is  significantly  different  from  that  of  sea  water 
sulphur  whose  6S34  =  20.4  °/00.   This  concept  is  strongly 
supported  by  the  work  of  Skerman  (1969)  who  demonstrated  co- 
populations  of  aerobes,  anaerobes  and  sulphate-reducers  in 
cores  taken  in  Milford  Sound,  New  Zealand.   He  found  these 
three  types  of  bacteria  co-existing  in  ratios  of  6:7:5  and 
6:6:5  in  the  upper  inch  or  two  of  the  sediments.   In  one 
such  core,  co-existent  populations  of  all  three  types  of 
bacteria  were  demonstrated  to  exist  down  to  7  ft  in  an  8  ft 
core.   The  depth  of  aerobic  activity  was  attributed  to  the 
strictly  terrestrial  origin  of  the  deposits.   Kaplan  and 
Rafter  (1964)  also  showed  that  surface  sediments  support 
Thiobacilli  who  draw  02  from  overlying  waters  and  H2S  from 
the  lower  sediments  by  diffusion. 

The  sulphates  in  the  Santa  Monica  and  Santa  Barbara 
cores  (Kaplan  e_t  a_l,  1963)  became  enriched  in  S32  with  depth 
while  the  Santa  Catalina  core  showed  a  progressive  enrich- 
ment in  S3I+  (Figure  16).   However,  this  data  was  scanty  and 
showed  only  a  trend  with  no  details. 

3 .   Isotope  Ratios  of  Sulphates  in  a  Core 

Kaplan  ejt  aJ  (1963)  suggested  a  hypothetical  model 
of  the  sulphur  cycle  in  recent  California  sediments.   They 
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Figure  15.   Sulphate  concentration  versus  depth  for  three 

cores  from  California  shelf  (after  Kaplan  et  al , 
1963)  . 


53 


100 


dS34S04    °/°o 
+10  +15  +20 


+25 


400-J 


Figure  16.   6S31*  values  for  three  cores  from  the  California 
shelf  (after  Kaplan  e_t  al ,  1963). 
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assumed : 

a.  Constant  rate  of  deposition. 

b.  Sulphate  was  present  in  the  sediment  column 
dissolved  in  the  interstitial  waters,  not  in  combined  form 
(barite,  gypsum)  or  adsorbed  onto  clays. 

c  .   There  was  no  diffusion  of  water  into  the 
sediments  from  the  overlying  water;  i.e.,  the  water  in  the 
column  is  a  closed  system  added  to  only  by  new  deposits  at 
the  top. 

d.  At  some  time  after  deposition,  oxygen  was  ex- 
hausted by  microbial  activity. 

e.  After  the  oxygen  was  all  depleted,  then  sulphate 
reduction  started  and  proceeded  to  completion. 

f.  A  cascade  of  reactions  then  occurred  to  sul- 
phides which  ultimately  gave  pyrite  as  an  end  product. 
The  consequences  of  this  model,  in  terms  of  what  happened 
to  the  sulphate  were  that: 

The  6S3lt  of  sulphate  sulphur  would  increase  with 
depth  and  exceed  that  of  sea  water  at  all  depths 
after  the  initiation  of  sulphate  reduction. 
To  this  model,  certain  modifications  were  added  to  explain 
their  findings: 

a.   At  the  sediment-water  interface  sulphate  reduc- 
tion (pyritization)  could  occur  in  microenvironments  such  as 
faunal  tests.   This  pyrite,  formed  from  the  infinite  sulphate 
pool  of  the  overlying  water,  was  depleted  in  S 3  *•  and  could 
later  serve  as  a  source  of  light  sulphur  for  the  formation 
of  light  sulphate  by  sulphur-oxidizing  reactions. 
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b.   The  isotopic  enrichment  in  S 3  **  was  found  to  be 
inversely  proportional  to  the  rate  of  reduction.   That  is, 
the  faster  the  rate  of  reduction,  the  less  the  fractionation 
between  the  parent  sulphate  and  the  product  pyrite. 

With  regards  to  the  oxygen  ratios  in  sulphates,  the 
following  processes  are  therefore  suggested  as  likely  to 
occur : 

a.  At  the  sediment-water  interface  both  aerobic 
and  anaerobic  bacteria  can  co-exist  and  can  reduce  and  oxi- 
dize sulphur  in  relatively  close  proximity  to  each  other. 
However,  since  the  Thiobacilli  have  been  shown  not  to  frac- 
tionate elemental  sulphur  (Kaplan  and  Rittenberg,  1964; 
Jones  and  Starkey,  1957)  in  sulphate  production  and  sul- 
phides only  weakly  (Kaplan  and  Rittenberg,  1962;  Kaplan  and 
Rafter,  1958),  it  is  assumed  that  their  fractionation  effect 
upon  oxygen  is  also  small  or  non-existent.   This  is  perhaps 

a  poor  assumption  but  there  is  no  other  data  bearing  directly 
on  the  problem. 

b.  As  sediment  is  buried  more  deeply  in  the  core, 
its  oxygen  supply  becomes  more  and  more  depleted  by  the 
action  of  Thiobacilli  and  the  possibility  for  completely  re- 
ducing microenvironments  to  occur  would  increase.   Thus  the 
sulphate  is  being  more  depleted  in  016  and  becoming  heavier. 

c.  Nitrophilic  bacteria  could  oxidize  free  sulphur 
to  form  sulphates  using  nitrates  as  an  oxygen  source.   This 
reaction  would  probably  only  occur  after  all  the  oxygen  was 
used  up  since  free  oxygen  has  a  lower  bonding  energy  than 
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nitrates,  or  may  occur  to  an  increasing  extent  in  microenviron 
ments.   This  sulphate  so  formed  would  be  light  in  sulphur 
relative  to  the  elemental  sulphur  which  is  itself  light  rela- 
tive to  the  ambient  sulphates.   This  is  expected  since  each 
process  from  pyritization  to  sulphur  oxidation  results  in  a 
fractionation  of  the  sulphur  with  the  lighter  S32  being  pref- 
erentially acted  upon.   It  is  assumed  that  the  sulphate  so 
formed  would  similarly  be  lighter  in  oxygen  than  the  nitrates 
from  which  it  was  formed.   There  does  not  appear  to  have  been 
any  work  done  on  the  oxygen  exchange  time  constant  for  ni- 
trates.  Longinelli  and  Craig  (1967)  remarked  that  the  ni- 
trate time  constant  was  presumed  to  be  long  but  gave  no 
further  data.   It  is  presumed  that  it  is  at  least  as  long  as 
that  for  sulphates. 

d.  The  sulphate  oxygen  is,  at  each  step,  continuing 
the  slow  equilibration  process  with  the  water.   After  all 
other  processes  are  subtracted  out,  a  residual  decline  in 
the  oxygen  ratio  approaching  that  of  the  ambient  interstitial 
water  should  be  noted. 

e.  The  effect  of  organically  bound  sulphur  on  the 
sulphate  system  is  ignored.   Kaplan  et  al_  (1963)  estimated 
that  organically-bound  sulphur  constituted  less  than  \%    of 
the  total  sulphur  in  the  sediments. 

The  sum  total  of  these  processes,  except  time 
equilibration  is  presented  as  an  idea-gram  in  Figure  17. 
This  figure  shows  only  a  qualitative  hypothetical  model  and 
does  not  take  into  account  any  other  processes.   From  Figure 
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Region  I     Sulphate  reducers  (  H20  +  C02) ;  02  + 
sulphate  producers  (  -H20) 

Region  II    Sulphate  reducers  (  H20+C02);  N03  + 
sulphate  producers  (  -H20) 

Region  III   Sulphate  reducers  [step  II  of  Harrison  and  Thode 
(1958)  dominant] 

Region  IV   Sulphate  reducers  [step  I  of  Harrison  and  Thode 
(1958)  dominant] 

Figure  17.   Idea-gram  of  nutrients  and  oxygen  ratios  versus 

depth.   (Depth  scale  is  arbitrary  and  concentra- 
tions of  nutrients  are  arbitrary) . 
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17  it  can  be  seen  that  the  sulphate  concentration  decreases 
with  depth  while  being  reduced  and  oxygen  is  being  used  up 
until  nitrophilic  bacteria  start  producing  more,  after 
which  it  monotonically  decreases  in  concentration.   It  was 
noted  qualitatively  during  sample  handling  that  the  amount 
of  precipitate  in  the  cores  did  decrease  by  a  factor  of 
about  5  from  top  to  bottom  of  the  cores. 
4 .   Effect  of  pH  on  Oxygen  Ratios 

If  no  processes  acted  on  the  water  or  the  sulphates, 
and  the  sulphates  were  initially  heavier  than  the  water  by 
9.5  °/oo  (Longinelli  and  Craig,  1967),  then  one  would  expect 
the  oxygen  ratios  to  approach  each  other  in  a  closed,  con- 
stant volume  system  (Figure  18).   The  S018H  Q  and  SO18™ 
should  asymptotically  approach  each  other  and  in  time  have 
a  common  del  value  between  9.5  and  0.0  per  mi  lie .   Teis  (1956) 
noted  that  the  rate  of  equilibration  is  pH  dependent.   In- 
creased C02  concentration  speeds  the  equilibration.   Lloyd 
(1967)  noted  this  pH  dependence  and  ran  several  of  his  sam- 
ples at  different  pH ' s  to  determine  its  effect  on  the  time 
rate  of  equilibration.   From  Longinelli's  graph,  one  can 
determine  that 


d  (log  tjj 


d(pH) 


T°C  =  constant 


=  0.35 


For  instance,  the  half  life  Lloyd  calculated  for  pH  7  was 
about  3  years  (at  200°C)  and  was  50  years  at  pH  9  (at  200°C)  . 
Thus,  if  the  pH  in  the  sediment  column  is  to  change  with 
depth  (time) ,  the  rate  of  exchange  will  change  correspond- 
ingly.  In  all  levels  of  the  previously  discussed  model, 
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Figure  19.   Idea-gram  of  601 8  versus  depth  (with  pH  varia 
tion  similar  to  that  found  by  Kaplan  et_  al_, 
1963) . 
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C02  is  produced  by  bacterial  action.   This  will  quickly 
equilibrate  with  the  interstitial  water  as  well  as  change  its 
pH.   Kaplan  e_t  al_  (1963)  noted  an  initial  decrease  in  pH  in 
the  first  30-100  cm  of  core  in  their  Santa  Barbara,  Santa 
Monica  and  Santa  Catalina  basin  cores.   Below  that  level 
the  pH  was  erratic  but  generally  increased  slowly  with  depth 
to  a  value  of  about  8.0  (from  7.4).   This  can  be  interpreted 
to  mean  that  the  C02  being  added  to  the  core  in  situ  below 
100  cm  by  sulphate  reducers  was  not  sufficient  to  prevent 
the  slow  increase  of  the  pH  of  the  interstitial  water  in  the 
core  from  other  processes.   The  surface  layers  may  have  had 
so  much  sulphate  that  the  C02  produced  could  significantly 
lower  the  pH  as  was  reflected  in  their  analysis.   A  general- 
ized idea-gram  of  6018cn   versus  depth  is  drawn  in  Figure  18. 
The  effect  of  the  expected  variation  of  pH  in  a  core  upon  the 
rate  of  oxygen  isotope  equilibration  is  to  speed  it  up  in  the 
surface  layers  and  to  slow  it  down  at  depth.   If  it  varies  by 
0.5  pH  units,  the  time  constant  could  be  approximately  halved 
The  effect  of  pH  upon  the  del  curves  is  idealized  in  Figure 
19.   In  reality  the  effect  of  pH  on  the  shape  of  the  two 
isotope  ratio  curves  may  be  very  small  over  30-100  cm  of  core 

In  the  surface  layers  and  down  to  the  bottom  of  the 
modeled  nitrophilic  zone,  water  is  being  used  up  by  both 
nitrophiles  and  aerobic  sulphate  producers.   Again  one  ex- 
pects the  lighter  water  to  enter  the  reaction  first.   And 
we  note  that  sulphate  reducing  bacteria  respire  water  as  a 
by  product.   Thus,  exclusive  of  the  time  effect,  the  SO18,*  q 

61 


will  be  a  combination  of  these  two  biological  reactions 
(Figure  20) .   Figure  20  assumes  that  sulphate  reduction 
(water  production)  is  proportional  to  the  sulphate  concen- 
tration and  that  sulphate  production  (water  loss)  is  inversely- 
proportional  to  the  sum  of  the  oxygen  concentration  and  the 
nitrate  concentration:  i.e., 


d[H201   =  Kx  [SO,2"] 

a 


d[H20]b  -  K2  [02  +  N(V] 

Thus, 

d[EH20]  =  Kj  [SO,2-]  +  K2  [02  +  N032"] 

The  sulphate  present  in  the  top  layers  of  both  cores 
was  noted  as  being  lighter  or  at  least  the  same  as  the  sul- 
phates in  the  overlying  waters  (Figure  6)  and  thus,  water 
produced  from  their  reduction  would  be  lighter  than  the  am- 
bient interstitial  water  giving  rise  to  an  initial  decrease 
in  isotope  ratio  in  the  interstitial  waters  with  depth.   As 
the  lighter  sulphates  were  used  up,  the  water  became  heavier 
and  heavier.   Down  to  30  inches  of  core  depth,  data  from  core 
B-l  followed  this  hypothesis,  below  this  (40  -  50  inches)  it 
became  dramatically  lighter  as  though  some  new  mechanism  were 
then  dominant.   The  ratios  for  sulphate  (Figure  8)  also 
shifted  down  and  there  was  a  decided  loss  of  O18  from  the 
water- sulphate  system. 
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Figure  20.   Idea-gram  of  018(H20)  and  sulphate  concentration 
versus  depth  as  affected  by  sulphate- reducing 
bacteria . 
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5 .   Physical  Processes  Which  Act  on  the  Water-Sulphate 

System 

Other  processes  acting  on  the  water- sulphate  system 
independent  of  anything  thus  far  discussed  cannot  be  exactly 
specified.   Garlick  and  Dymond  (1970)  showed  that  volcanic 
glass  shards  do  exchange  O18  with  water  in  sediments.   They 
showed  that  probably  the  glass  exchanges  oxygen  very  slowly, 
altering  the  del  value  +  12  °/0o  in  49  million  years.   Savin 
and  Epstein  (1970b)  showed  that  there  was  no  evidence  of  ox- 
ygen isotope  exchange  in  detrital  clay  minerals  from  24 
worldwide  ocean  cores.   However,  authigenic  montmorillonite 
was  shown  to  fractionate  oxygen  relative  to  water  in  the 
ratio  of  1.028  to  1.000.   They  stated,  "when  the  minerals 
of  an  ocean  sediment  are  formed  in  situ,  this  may  be  de- 
tected by  their  isotopic  composition"  (being  different  from 
the  same  mineral  of  detrital  origin) .   They  showed  that 
minerals  formed  in  a  marine  environment,  "suggest  formation 
in  the  marine  environment  under  isotopic  equilibrium  con- 
ditions."  In  point  of  fact  mixtures  of  minerals  can  be  both 
lighter  or  heavier  isotopically  then  their  detrital  counter- 
parts.  Siever  (1968)  showed  that  different  clays  recon- 
stituted in  silicate-spiked  water  could  become  less  hydrous. 
Boron  also  has  been  shown  to  fractionate  on  adsorption  to 
clays  as  boric  acid  (Schwarz,  Aguei  and  McMullen,  1969). 

Kaplan,  Sweeney  and  Nissenbaum  (1969)  analyzed 
several  cores  from  the  Red  Sea  hot  brine   deeps  for  sulphur 
compounds  and  sulphur  ratios.   Figure  21  is  a  histogram  of 
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Figure  21.   Histogram  of  sulphur  isotope  distribution  of 
different  sulphur  species  from  interstitial 
water  and  sediment  from  Atlantis  II  Deep  in  the 
Red  Sea  (cores  84k  and  120k)  (after  Kaplan, 
Sweeney  and  Nissenbaum,  1969). 
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their  data  from  those  cores  taken  from  the  Atlantis  II  Deep 
whose  sediments  have  been  shown  to  be  sterile  (Watson  and 
Waterbury,  1969;  Triiper,  1969).   The  interstitial  sulphate 
was  being  fractionated  when  it  entered  the  solid  phase  and, 
further,  the  sedimentary  and  interstitial  sulphates  were 
being  fractionated  by  the  physical  processes  which  change 
them  to  sulphides  and  elemental  sulphur.   These  processes 
were  apparently  step  processes  as  a  range  of  values  for  S34 
were  not  found  in  the  cores.   In  any  event  some  magnitude 
of  the  changes  to  the  sulphur  population  can  be  obtained 
from  this.   It  is  expected  that  the  oxygen  population  would 
show  corresponding  responses  to  physical  processes. 

The  effects  of  all  of  the  processes  discussed  is 
diagrammatically  depicted  by  Figure  22.   There  are  three 
processes  which  are  capable  of  adding  016  to  sulphates  in 
the  interstitial  water  of  a  core  section;  they  are: 

a.  Biological  oxidation  of  H2S  to  sulphate  by 
aerobic  bacteria 

b.  Biological  oxidation  of  H2S  to  sulphate  by 
nitrophilic  bacteria 

c.  Time  dependent  equilibration  with  water 


H2016  +  S0316018   ->  HoO18  +  SO, 


1  6 


There  are  two  processes  which  can  remove  016  from 
the  sulphates  of  a  core  sample;  they  are: 

a.  Biological  reduction  of  sulphates 

b.  Chemical  oxidation  of  H2S  (aerobically) . 
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Figure  22.   The  water- sulphate  oxygen  isotope  system  in 
sediments.   (Arrows  indicate  preferential 
movement  of  01 6 . ) 


67 


These  same  five  reactions  operate  on  water  in  an  inverse 
fashion  to  alter  the  oxygen  isotope  ratio  Ru  n.   In  addition, 

112'-' 

other  processes  can  operate  to  add  or  subtract  016  from  the 
water.   Examples  of  these  were  discussed  above. 


B.    DEVELOPMENT  OF  MATHEMATICAL  MODEL 

The  possible  effect  of  nitrophilic  bacteria,  diffusion, 
reduction  of  organic  sulphur  compounds  and  very  low  sulphate 
concentrations  will  not  be  considered  further.   It  was  as- 
sumed that  in  core  B-l  the  zone  10  to  50  inches  was  anaerobic; 
thus,  aerobic  reactions  were  excluded.   This  left  biological 
reduction  of  sulphates,  time  equilibration  of  water- sulphates 
and  other  processes  to  be  considered  in  developing  a  mathe- 
matical model. 

1 .   Effect  of  Sulphate-Reducing  Bacteria 

Kaplan  and  Rittenberg  (1964)  showed  that  the  amount 
of  fractionation  of  sulphate  S31*  is  inversely  proportional 
to  the  rate  of  reduction  and  independent  of  sulphate  con- 
centration.  Such  a  relationship  was  presumed  to  hold  for  ox- 
ygen ratios : 


3R 


SOi 


3t 


3R 


biol   "  klR 


H20 


SO, 


3t 


biol 


(1) 


where 


ki  =  proportionality  constant. 


2 .   Effect  of  Other  Processes 

The  rate  of  change  of  the  oxygen  ratio  in  water  due 
to  other  processes  was  arbitrarily  assumed  to  be  proportional 
to  the  rate  of  clay  formation  and  had  a  fractionation  factor 
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of  1.028  (after  Savin  and  Epstein,  1970a) 


3R 


H,0 


3t 


o.p 


x1  3c 
=  ±k2-^  w 


(2) 


where 


c  =  concentration  of  authigenic  materials 

w  =  dx/dt  =  rate  of  compaction  rate  of  deposition 
k2  =  1.028. 

3 .   Effect  of  Time  Equilibration 

The  rate  of  equilibration  of  sulphate  with  water 

was  taken  to  be  proportional  to  the  difference  of  the  oxygen 

ratios  of  water  and  sulphate: 
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where    k3  =  time  constant. 


4 .   Computational  Model 

Now  some  expression  was  needed  to  express  the  con- 
centration of  authigenic  minerals  with  depth.   It  was  assumed 
that  the  change  of  concentration  with  time  was  proportional 
to  the  concentration  of  free  interstitial  water  present, 
as  suggested  by  the  necessity  of  the  presence  of  water  for 
authigenic  mineral  formation  (Siever,  1968;  Savin  and  Epstein 
1970a).   Further,  the  rate  of  compaction  plus  the  rate  of 
deposition  was  assumed  constant  with  depth: 


dx 

dt 


constant 


3c   dx  _  3c 

at  dt  "  at 
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and 


It  ■  k^  tH2o]  (4) 


at 

where  k^  =  proportionality  constant. 

Figure  23  represents  the  moisture  content  of  a  core  taken 
by  Kaplan  et  al  (1963) .   It  can  be  seen  that  the  moisture 
content  decreases  with  depth.   The  Santa  Barbara  core  (7131) 
had  a  nearly  constant  composition  of  clay/silt  with  depth 
and  thus  probably  represented  as  uniform  a  core  as  one  could 
obtain.   The  clay/silt  content  of  core  A-2  (Figure  24)  was 
nearly  the  same  as  that  of  the  Santa  Barbara  core. 

An  exponential  curve  was  fitted  to  the  Kaplan  data 
presuming  that  the  water  content  was  being  progressively 
used  up.   The  equation  for  this  curve  is: 


where 


ru  nl     -0.00774x    , 
[H20]  =  e  gm/gm 


x   =  depth  in  cm 

K   =  0.774xl0"2  . 
w 


3c   ,    -Kw*x 
7TT-  =  kit  e 


Referring  to  equation  (4)  one  can  write 

i£ 

3t 

and  substitution  into  equation  (2)  gave 

— fe^       =  k2  k,  e   w  (5) 

dt    O.p. 

Combining  (1)  and  (3)  to  give  the  total  time  rate  of  change 
of  the  sulphate  ratio  gave 

■fe^  ■ k-  dr  -  k3  (Rso.  -  rh2o^  w 


70 


WATER     CONTENT 
(%   of  dry  weight) 
100  200  300 


Figure  23.   Moisture  content  of  a  core  (after  Kaplan  et  al, 
1963)  . 
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Figure  24.   Clay  versus  silt  for  core  A-2 
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Similarly,  for  water  ratios,  combining  equations  (1) , 
(3)  and  (5)  gave 

dRH   0  1  ~VX 

— at-  =  "  kl  RZT  +  k2  k*  e  "  k3   ^RS0,  "  RH2(P 

(7) 

Equations  (6)  and  (7)  are  expressions  representing 
the  effects  upon  the  water- sulphate  system  by  bacteria,  phy- 
sical processes  (authigenic  mineral  formation) ,  and  time  rate 
of  equilibration.   These  equations,  when  added  together,  give 
a  computational  equation  to  derive  the  time  constant  sought, 

dRqp.     ^^u  n  "   *x 

-jSOl  +  _H|£  .  fc2  ^    e  «    .  2k3  (R^  .  Rh2q)   (8) 

The  analysis  of  core  A-2b  showed  that  there  was 
little  likelihood  that  a  significant  portion  of  the  018  in 
the  water-sulphate  system  was  being  lost  to  the  formation  of 
authigenic  minerals  or  other  physical  processes.   Therefore, 
the  first  term  on  the  right  side  of  equation  (8)  was  dropped 
leaving : 

dR„n     dR„  n 

in  which  k3  is  the  looked  for  time  constant  for  sulphates  in 

-  i 
years 

Berner  (1964)  assumed  that  the  deposition  rate  for 

recent  California  shelf  sediments  was  about  0.30  cm  per  year. 

If  this  rate  is  assumed  to  be  approximately  true,  then  each 

10  inches  of  core  segment  represented  84.5  years.   This  was 
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adopted  as  the  time  scale  for  core  B-l.  A  pelecypod  shell 
was  found  at  the  botton  (50  inches)  of  core  B-l  and,  if  it 
were  possible  to  have  it  analyzed  for  age  by  Clh  analysis, 
the  time  scale  assumed  above  could  be  revised  accordingly. 
However,  such  an  analysis  was  not  available.  Core  B-l  then 
was  assumed  to  represent  a  time  series  422.5  years  old. 
This  was  incorporated  in  Figure  25. 

With  each  time  step  being  84.5  years,  the  computa- 


tional model  became 
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H20 
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SO, 
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H20 


+RH20 


=  k3  years 


-  i 


(10) 


Equation  (10)  utilized  simple  difference  molecules  for  the 
time  rate  of  change  of  the  sulphate  and  water  ratios  and 
arithmetic  averages  for  the  absolute  ratio  differences.   The 
computed  k3's  are  presented  in  Table  9.   The  first  10  inches 
of  the  core  was  apparently  undergoing  processes  which  were, 
as  expected,  contrary  to  the  assumptions  made  in  the  mathe- 
matical model  and  thus  was  disregarded.   The  average  of  the 
rest  of  the  data  gave  a  value  of 


k3=  -0.0150±0.0029(standard  deviation) years 


-  i 


The  value  of  k3  is  negative  since  the  ratios  decrease  each 
year . 
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Figure    25.      Combined  parameters    of    core    B-l 
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Table  9 
Values  of  Time  Constant  Derived  From  Evaluation  of 

Equation  (9) 


Core  Depths 
(Inches) 
(Core  B-l) 

0  -  2 


10 
20 
30 
40 
50 


RSCU    RH20    Time  Constant  (k3) 


42.2 
43.4 
43.5 
43.5 
43.0 
42.7 


24.1 
23.6 
25.4 
26.5 
23.0 
21.7 


0.00465 
0.0134 
0.0118 
0.0151 
0.0195 


Average  0  .  0150±0 . 0029  years"1 
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5  .   Calculation  of  t-g 

Referring  back  to  equation  (3),  the  value  of  the 
half  life  (th)    of  the  sulphate  oxygen  was  finally  computed, 


dR 


SO 


dt 


time  "  k3  (RS0,  "  RH20} 


The  solution  to  this  ordinary  differential  equation  is 
ln  (RS04  "  RH2(P  =  k3  At 


(3) 


(ID 


When  the  argument  of  the  log  term  on  the  left  is  h,    then 
t=th   and, 

th   =  46  .  2  years 
This  is  approximately  two-thirds  the  value  which  Teis  (1956) 
obtained  of  70  years  and  three  orders  of  magnitude  less  than 
what  Lloyd  (1967)  obtained.   The  figure  of  46.2  years  is  in- 
versely proportional  to  the  estimate  of  the  time  scale  assumed, 
Therefore,  it  appeared  that  Teis'  estimate  was  closer  to  what 
may  be  found  in  the  real  ocean. 

The  Lake  Vanda  study  (Rafter  and  Mizutani,  1967) 
showed  that  the  sulphates  were  changing  their  oxygen  isotope 
ratios  more  rapidly  than  was  the  water.   This  may  indicate 
that  the  physical  processes  acting  on  the  water  were  opera- 
ting in  the  opposite  direction  then  the  other  biological 
processes . 

Further  work  remains  to  define  more  accurately  the 
oxygen  cycle  in  sulphates  in  sediments  and  the  ocean  so  that 
more  accurate  relationships  than  equations  (6)  and  (7)  can 
be  formulated. 
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V.   SUMMARY 

The  original  intent  of  this  paper  was  to  determine  the 
time  constant  for  the  rate  of  equilibration  of  oxygen  in 
sulphates  with  water.   This  was  done  experimentally  using 
two  cores  from  the  upper  Carmel  Canyon  as  a  time  series. 
The  interstitial  sulphates  and  the  interstitial  waters  were 
analyzed  for  their  oxygen  isotope  ratios  on  a  mass  spectro- 
meter.  Next,  the  various  biological  and  physical  processes 
which  can  act  on  sediments  to  alter  the  water- sulphate  sys- 
tem isotopically  were  reviewed  and  discussed.   These  pro- 
cesses were: 

a.  Biological  reduction  of  sulphate  by  sulphate- 
reducing  bacteria  (Desulf  ovibrio  sj) . )  . 

b.  Biological  oxidation  of  sulphur  compounds  by 
aerobic  sulphate-reducing  bacteria  (Thiobacilli  sp . ) . 

c.  Biological  oxidation  of  sulphur  compounds  by  an 
aerobic  sulphate  producing  bacteria  (e.g.,  Thiobacilli 
denitrif icans) . 

d.  Chemical  oxidation  of  sulphur  compounds  aerobically 

e.  Other  physical  processes  which  act  to  alter  oxygen 
isotope  ratios  of  water  which  do  not  affect  the  sulphate 
oxygen  ratios  (e.g.,  authigenic  mineral  formation). 

f.  The  time  rate  of  equilibration  of  the  oxygen  ratios 
of  the  water-sulphate  system. 
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A  mathematical  model  was  proposed  which  had  terms  modeling 
those  processes  considered  significant  in  the  anaerobic  zone 
of  a  sediment  column.   These  processes  were  a,  e  and  g  above. 
The  two  differential  equations  thus  obtained  were  manipulated 
to  yield  a  computational  model  which  gave  the  time  constant 
for  the  oxygen  ratios  of  sulphate  of, 

-0 .0150+0. 0029(standard  deviation)  years"1 

This  was  then  used  to  determine  the  half  life  of  sulphate  ox- 
ygen ratios  as 

th  =  46.2  years 
which  was  in  fair  agreement  with  Teis'  (1956)  work  but  three 
orders  of  magnitude  different  from  Lloyd's  (1967)  work.  Once 
such  constants  can  be  precisely  determined  as  a  function  not 
only  of  time  but  of  pH  and  biological  processes  as  well,  in- 
vestigators will  have  another  powerful  dating  tool  available 
for  water  mass  age  determinations  which  will  not  require  the 
use  of  radioactive  nuclide  techniques. 
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